Unified picture for diluted magnetic semiconductors 
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For already a decade the field of diluted magnetic semiconductors (DMS) is one of the hottest. 
In spite of the great success of material specific Density Functional Theory (DFT) to provide ac- 
curately critical Curie temperatures (To) in various III-V based materials, the ultimate search for 
a unifying model/theory was still an open issue. Many crucial questions were still without answer, 
as for example: Why, after one decade, does GaMnAs still exhibit the highest To? Is there any 
intrinsic limitations or any hope to reach room temperature? How to explain in a unique theory 
the proximity of GaMnAs to the metal-insulator transition, and the change from RKKY couplings 
in II- VI materials to the double exchange regime in GaMnN? The aim of the present work is to 
provide this missing theory. We will show that the key parameter is the position of the Mn level 
acceptor and that GaMnAs has the highest Tc among III-V DMS. Our theory (i) provides an overall 
understanding, (ii) is quantitatively consistent with existing DFT based studies, (iii) able to explain 
both transport and magnetic properties in a broad variety of DMS and (iv) reproduces the Tc 
obtained from first principle studies for many materials including both GaMnN and GaMnAs. The 
model also reproduces accurately recent experimental data of the optical conductivity of GaMnAs 
and predicts those of other materials. 



In recent years the so called diluted magnetic semi- 
conductors (DMS) have attracted much interest owing 
to their potential for spintronic devices [ij, |2| . As far as 
the symbolic barrier of room temperature can be crossed, 
their interest is to combine both the traditional electronic 
functionality and the spin degree of freedom. DMS are 
obtained by doping a semiconductor host with a little 
amount of transition metal (TM) ions. Ferromagnetism 
is mediated by carriers antiferromagnetically coupled to 
localized spins. Transport and magnetic in DMS proper- 
ties are very sensitive to dilution effects and to the pres- 
ence of intrinsic defects. Available theoretical studies 
can be divided into two different families (i) first prin- 
ciple (DFT) and (ii) model. Note that there are still 
many dissensions among these numerous studies. First 
principle based studies were successful to reproduce ac- 
curately Curie temperatures for both as grown and an- 
nealed samples Although powerful, first principle 
methods are essentially fully material-dependent. Thus 
it is difficult to draw very general conclusions and identify 
the most relevant physical parameters that control both 
magnetic and transport properties. Besides DFT studies, 
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the most used model is based on a six or eight realistic 
bands Kohn-Luttinger (KL) Hamiltonian [J includ- 
ing a pd-exchange interaction between localized spins and 
itinerant holes. However, (i) the pd-coupling is treated 
perturbatively and (ii) the dilution effects are neglected: 
this describes the valence band scenario. These studies 
are inconsistent with first principle calculations which, 
for example, show the existence of a unambiguous pre- 
formed impurity band (IB) in GaMnAs even at relatively 
high Mn concentration (4-7%). In addition the VB pic- 
ture for GaMnAs is unable to explain recent optical con- 
ductivity measurements in both as grown and annealed 
samples • Finally, those approaches are restricted to 
metallic systems and can not explain the variations of the 
Curie temperatures with both hole and magnetic impu- 
rity concentrations in III-V compounds. The aim of the 
manuscript is to propose a model which is able to pro- 
vide a overall understanding of the magnetic and trans- 
port properties both qualitatively and quantitatively in 
the whole family of III-V materials. The model Hamil- 
tonian (V-J model) that describes the carriers (holes or 
electrons) interacting with the localized impurity spins 
and it reads, 

H = - £ *«4rc,v + Yl JA - s * + J2 v A a a° (!) 

i icr 

The hopping term ty=t for nearest neighbors only. cj CT 
(resp. Cia) is the creation (annihilation) operator of a 



hole of spin a at site i. is the local coupling between 
localized impurity spin Si (S=5/2 for Mn 2+ ) and a spin 
carrier Sj(p band). Whereas the on-site potential V$ de- 
scribes the effect of the substitutional disorder. J^=Pi 3 p d 
and Vj =PiV where p^=l for impurity sites, otherwise 
pi=0. The magnetic impurities are randomly distributed 
on a simple cubic lattice and x denotes the impurity con- 
centration. Note that this minimal one band model was 
already shown to be a good candidate to describe qual- 
itatively the magnetic properties of various III-V DMSs 
0? ]. 

We now discuss how we fix this set of 3 parameters. 
In III-V hosts the bandwidth is rather close to that of 
GaAs (W^7-8 eV), therefore, for simplicity, we assume 
the same value t« 0.7 eV for all considered (III,Mn)-V 
compounds. Note that a variation of about ±10% of t will 
not affect considerably our results. In both (II,Mn)VI or 
(III,Mn)-V DMS J pd is about 1 eV [13, El, thus in the 
following we set J p( z=l-2 eV (accepted value for GaM- 
nAs). Thus, a single parameter is left (V). This param- 
eter V is set to reproduce for each compound the bound 
hybridized pd states energy (E&) with respect to the top 
of the VB. Note that this bound-state is three-fold de- 
generate, this is illustrated in the left panel of FigfSJ At 
finite, but small concentration on Mn and because each 
Mn provides a single hole this will lead to an impurity 
band 1/3 hole filled (see right panel of Figi2|) . In 

the case of a single Mn in GaAs the bound states lies 
at E h 110 meV 0, EH above the VB. In FigfTJ the 



calculated Et is plotted as a function of V/t. In this fig- 
ure the realistic values of Mn acceptor level in various 
hosts (GaAs,InAs,InN..) are shown. These values were 
obtained either experimentally or within first principle 
calculations [l7-21[. For GaMnAs we have found V=1.8 
t and 5 t for GaMnN. The inset of FigfTJ shows the cal- 
culated spin splitting A(x) = E^^ — E^ ax for GaMnAs 
as a function of x, E° nax is the largest eigenvalue in the 
a sector (see right panel of Fig [2] for an illustration). We 
see that our results are in excellent agreement with those 
obtained from ab initio LSDA calculations [22|, [23[ (open 
symbols). The dashed line shows the Mean Field Vir- 
tual Crystal Approximation (MF-VCA) value A=x3 p dS. 
In contrast to what is often found in the literature, the 
MF-VCA expression can not be used to extract J p( j in 
III-V DMS. For that reason the value that was found 
in [22j was much larger (4 times) than the value of 1.2 
eV. However, A=xJ p dS is a good approximation in II- 
VI materials as ZnMnTe or CdMnTe (small \V\ values). 
In our one band model each Mnc brings a single state. 
Thus, to be consistent with the realistic compound in 
which each Mn provides a single hole and n; = 3 p-d 
states, our model calculations for well annealed samples 
will performed at the hole density p — x/3. We will get 
back to this important point in the following. 

Let us now first discuss the procedure to calculate the 
magnetic properties in any compound characterized by 
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FIG. 1: (Color online) Calculated Mn acceptor level Eb (eV) 
as a function of V/t in III-V hosts (assuming t~ 0.7 eV and 
J P dS=4.3 t). The realistic Ef, [a-e] are extracted from [l7l - [2l| ]. 
(Inset) Spin splitting A (eV) as a function of the Mn con- 
centration x for Gai-xMn^As. Symbols correspond to LSDA 
calculations ( [i] [22] , [ii] [13| ) and the continuous line to the V- J 
model. The dashed line is the MF-VCA expression. 
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FIG. 2: (Color online) This cartoon only shows the ni — 3 
pd-states near the top of the valence band for well annealed 
GaMnAs sample (see left panel). In GaMnAs, each Mn 2+ 
brings a localized spin (S=5/2) and one hole, thus at finite 
concentration x the hole density is p = x. To be consistent 
with the 1/3 filled impurity band, within our one band model 
the calculations are performed for the hole density p — x/3 



its set of 3 parameters. For each configuration of disorder 
(distribution of magnetic impurities in the host) we diag- 
onalize exactly the Hamiltonian (TTJ in both spin sectors 
assuming fully polarized d-spins (for details see ref . [jjj ) . 
Typically we have used simple cubic systems size from 
16 3 to 24 3 sites and the average are done with a few hun- 
dreds of disorder configurations. The diagonalization for 
each disorder configuration (denoted c) provides the full 
spectrum, eigenvalues and eigenvectors denoted {E£ a , 
l*)<r a) ( a denotes the eigenstate index) needed to cal- 
culate magnetic couplings and transport properties. The 
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FIG. 3: (Color online) Curie temperature (in K) for 
Gai-zMn^As and Gai-^Mn^N as a function of Mn concen- 
tration x within both ab-initio [E| (open symbols) and model 
calculations (filled symbols). 



magnetic couplings between two localized spins, respec- 
tively at site i and j, is given by the generalized suscep- 
tibility 
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where the Green's function are 



:\3<t). 



Within our 



Tr^Gy^Gjj^du: 



(2) 

G IM = 
the local ex- 

As mentioned 



model 



uj — H+ie 

change splitting reduces to Ej 
before, since our model provides a single state per 
impurity (instead of 3 p-d states in the realistic ap- 
proach) then the magnetic exchange that have to be 
used to calculate the Curie temperature are defined 
as Ji.j{x,p) = niJij(x,p — p/ni). We note that the 
same argument applies for the transport properties. In 
other words, the optical conductivity will be defined as 
o{lo,p) = nia(w,p/ni) where a(u>,p/ni) is calculated 
within our one band model, the results will be discussed 
in the last section. Note that the overall shape of the 
couplings agree well with those obtained from first 
principle studies. For example, for Gai-^Mn^As they 
are rather short range and essentially ferromagnetic for 
well annealed compounds 24, 25] . To calculate the Curie 
temperature, the effective dilute Heisenberg Hamiltonian 
Hueis = —J2ijJij( x iP)Si • Sj (sums runs over sites 
occupied by Mn 2+ ) is solved within the self consistent 
local RPA theory (SCLRPA)^. Note also that we 
could also calculate other magnetic properties such as 
magnetization as a function of T, magnons spectrum, 
spin stiffness ... but this is not the scope of the present 
work. Note also that the accuracy and reliability of 
SC-LRPA was already demonstrated several times. 

In Fig|3] we have plotted the Curie temperature (in 
Kelvin) for both Gai_ x Mn x As and Gai-^Mn^N as a 
function of Mn concentration assuming well annealed 



i i i 1 1 1 1 i 

O V-J model 

□ [a] Olb] 
AM Old] <W 





E b (eV) 



FIG. 4: (Color online) Curie temperature (in K) as a func- 
tion of the bound state Eb (eV) and for a Mn concentration 
x = 5% in well annealed samples. The filled symbols cor- 
respond to our model, the open symbols to Tc's computed 
with ab initio exchange integrals ([a-d] are respectively from 
Refs.&Sll.IIl). 



samples (p=x). Surprisingly, the agreement between 
model and first principle based calculations [U, @, [26| is 
not only qualitative but also quantitative in the over- 
all range of Mn concentration. It is important to stress 
again that no fitting parameter is used. Note that the 
ab initio results were themselves in very good agreement 
with the experimental data 27 - 29[ . These results already 
clearly demonstrate that our model captures accurately 
the physics in these materials, this will be confirmed in 
the following. Note that the nature of the couplings 
are very different in these two compounds. Indeed, in 
GaMnN the couplings are of double exchange type whilst 
more extended in the case of GaMnAs. Note also that 
a variation of about ±10% of the hopping t does not af- 
fect considerably the quantitative results (box symbol for 
x = 0.05). 

In FigH] we now show the model calculations of the 
Curie temperatures (in K) as a function of Et (eV) for 
various (III,Mn)V compounds. We present only the case 
of x = 5% of substitutional Mn and a hole density p=x. 
In this figure, available Tp's obtained from ab-initio ex- 
change integrals 0, [26|, 30, 31 1 are also plotted. Again 
we obtain a very good agreement between the V-J model 
calculations and those obtained starting from first prin- 
ciple approaches. The model calculations show that the 
Tc(Efc) curve exhibits a clear and pronounced resonant 
peak structure and predicts that Gai-^Mn^As has the 
highest critical temperature among III-V Mn doped ma- 
terials. By mean of transport calculations, the mobility 
edge is found close to the value V=2.1 t (E&w 0.15 eV). 
Interestingly, the narrow resonant peak (30 < Eh < 200 
meV) appears to be located near the metal-insulator 
phase transition. Thus, this already indicates why exper- 
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imentally both metallic and insulating compounds were 
obtained and why after annealing as-grown samples the 
materials exhibit an insulator-metal transition. For small 
E b (or V), the couplings are RKKY like and thus the 
Curie temperatures are very small 24j, |32j . This region 
corresponds to II- VI materials as ZnMnTe or CdMnTc. 
As Eb increases the critical Curie temperature gets larger 
resulting from the suppression of the RKKY oscillations 
and the couplings become shorter range (resonant effects 
due to the preformed impurity band). By increasing fur- 
ther V (thus Ef, ) the range of the couplings becomes 
shorter and shorter and the relevant couplings that con- 
trols the critical temperature get smaller, thus Tc re- 
duces too. For very large E}, the couplings become dou- 
ble exchange like leading to much smaller Tc. In this 
region the physics of percolation becomes more crucial 
(included in our theory). This is the case of GaMnN 
or AlMnN in which the nearest neighbor coupling dom- 
inates strongly. Within our theory Gao.95Mno.05P is 
found to be insulator in agreement with experimental 
observations [jlj . Note also that our predictions for 
Gao.95Mno.05P and Ino.g5Mno.05As (respectively ~ 50 K 
and 60 K) are close to the experimental maximum mea- 
sured value (T^ xp £=; 60 K for annealed samples with 

x « 6%) jiiilaJE^. 

In this section we now discuss whether our theory 
is also successful to explain transport measurements. 
We have performed the optical conductivity as a func- 
tion of the frequency u> in well annealed Ini-^Mn^As, 
Gai-^Mn^As and Gai_ x Mn x P for x = 5%. As discussed 
above, the optical conductivity is defined by a(u>,p) = 
n i So- ^a-( u} ;P/ n l)- The regular part of the one band op- 
tical conductivity is given within the Kubo formalism by, 



W a (w, p/m) = - — (n£ - nf)AZ tfl 5(hu ~E° a +EJ 3 ) 

a 

(3) 

N = L 3 is the total number of sites, a is the simple cu- 
bic lattice parameter. The lattice parameter of GaAs fee 
lattice is ao ~ 5.65 10 -10 m (4 atoms per unit cell), to 
keep the same unit cell volume a is state to a = -ffis- 
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where j x is the 
current operator in the x direction and is the occupa- 
tion number of the state \^%). In FigJ5]we compare the 
results of our model calculations with available experi- 
mental data for Gao.95Mno.05As [7]. We find an excellent 
quantitative agreement on the whole frequency range. In 
particular the peak location at 0.2 eV is accurately repro- 
duced. Note that we have recently performed a detailed 
analysis of the effects of compensating defects on the 
transport properties of GaMnAs [3(| in presence of com- 
pensating defects and it leads to excellent quantitative 
agreement with recent experimental data 0] ■ In particu- 
lar our theory is in agreement with the observed red-shift 
of the broad peak while increasing the hole density from 
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FIG. 5: (Color online) Model calculation of the conductivity 
(in unit fi" 1 ■ cm" 1 ) as a function of uj/t for well annealed 
(p = x) Ino.95Mno.05As, Gao.95Mno.05As and Gao.95Mno.05P- 
The symbols correspond to experimental measurements in 
Gao.95Mno.05As 0. 



compensating to annealed GaMnAs samples. This fea- 
ture is in agreement with the fact that Gao.95Mno.05As 
exhibits a preformed impurity band. In the same fig- 
ure, we also predict some trends for Gao.95Mno.05P and 
Ino.95Mno.05As. At low frequency the conductivity of In- 
MnAs is seen to be higher than the two other compounds 
which is correlated with the fact that the Fermi level in 
InMnAs is deeper in the valence band (more extended 
states), let us recall that the Mn level acceptor lies at 
Efc « 300 meV. For this material the optical conductiv- 
ity peak is located at 0.11 eV. For Gao.95Mno.05P the 
shape of <t{uj) exhibits two clearly distinct peaks. The 
principal one is located at u> ~ 0.2 eV and the narrower 
secondary peak lies at ui ~ 1.2 eV. For this Mn density it 
is interesting to notice that GaMnP is at the edge where 
the impurity band separates from the valence band. It's 
secondary peak originates from electronic transitions be- 
tween the separated impurity band and the top of the 
VB. As we move to larger Ef, we expect a shift of this 
secondary peak to highest frequency (case of GaMnN for 
example). For these last compounds it should be interest- 
ing to have transport measurements in order to confirm 
or infirm these predictions. 

To conclude, we have drawn up a model for diluted 
magnetic semiconductors based on the minimal V-J 
Hamiltonian treated in a non perturbative way and with- 
out the use of effective medium. The parameters of the 
model are fixed in agreement with first principle band cal- 
culation. For example, concerning well annealed GaM- 
nAs (and without inhomogeneities), that process allowed 
us to reproduce quantitatively the Curie temperatures 
obtained from first principle study which were found to 
be in excellent agreement with the experimental values. 



■5 



While varying the local potential V (or equivalently the 
bound state E^), the resonant structure of the critical 
temperature shows that GaMnAs can be considered as 
optimal among the (III,Mn)-V family. Moreover our the- 
oretical model is in excellent agreement with recent mea- 
surements of the optical conductivity for GaMnAs. For 
the band gap DMS GaMnN, the critical temperatures ob- 
tained from our model are also in good agreement with 
first principle calculations. For the first time a theoretical 
model allow to describe qualitatively and quantitatively 
both transport and magnetic properties of a broad range 
of diluted magnetic semiconductors and on both sides of 
the metal-insulator phase diagram. This unifying pic- 
ture bridges the gap between model approach and first 
principle calculations. Successes of this model, show that 
it represents a tool of choice to understand experimen- 
tal data. In addition it allows interesting predictions for 
specific materials. 
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